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ABSTRACT 

O 

Cool weakly ionized gaseous rotating disk, are considered by many models as 
the origin of the evolution of protoplanetary clouds. Instabilities against pertur- 
bations in such disks play an important role in the theory of the formation of 
stars and planets. Thus, a hierarchy of successive fragmentations into smaller and 
smaller pieces as a part of the Kant-Laplace theory of formation of the planetary 
system remains valid also for contemporary cosmogony. Traditionally, axisym- 
metric magnetohydrodynamic (MHD), and recently Hall-MHD instabilities have 
been thoroughly studied as providers of an efficient mechanism for radial transfer 
of angular momentum, and of density radial stratification. In the current work, 
the Hall instability against nonaxisymmetric perturbations in compressible ro- 
tating fluid in external magnetic field is proposed as a viable mechanism for the 
azimuthal fragmentation of the protoplanetary disk and thus perhaps initiating 
the road to planet formation. The Hall instability is excited due to the combined 
effect of the radial stratification of the disk and the Hall electric field, and its 
groth rate is of the order of the rotation period. Such family of instabilities are 
introduced here for the first time in astrophysical context. 
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Introduction 



Linear mode analysis provides a useful tool for gaining important insight into the rel- 
evant physical processes that determine the stability of rotating fluid configurations. The 
importance of magnetic fields in rotating disks has been demonstrated by the rediscovery 
of the magneto rotational instability (MRI) in which hydrodynamically stable flows with 
an angular velocity t hat is decreasing outwa rds are highly unstable when threaded by a 



weak magnetic field (IBalbus fc Hawley 19911 ). That investigation has been carried out in 



the magnetohydrodynamic (MHD) limit and invoked a number of approximations appro- 
priate to the study of the evolution of long-wavelength perturbations in the weak magnetic 
field limit. The inclusion of the Hall electric field (Hall MHD) is a relatively recent devel- 



opment (IWardle 1999: 


Ba 


bus & Tero 


uem 2001: 


Sano & Stone 2002a: 


Sano & Stone 2002bl: 


Salmeron & Wardle 2003: 


Desch 2004 


:llJrpin & Rudieer 2005: 


Rudieer & Kitchatinov 20051) 



The Hall electric field plays an important role in the disk's dynamics when the coupling be- 
tween the electrons and neutral components of the fluid is low. In such cases, the inertial 
length of the ions is longer than the characteristic perturbation's length scale and conse- 
quently the motions of the ions and electrons are decoupled. Indeed, it has been shown that 
the Hall electric field has a profound effect on the structure and grouth rate of unstable 
modes like the MRI's. Furthermore, as will be shown below, the Hall term gives rise to 
new branches of unstable modes. In particular, it will be shown that the Hall term, in the 
presence of radial stratification, excites non axisymmetric instabilities. 

As an astrophysical interest, we mention magnetically supported cool molecular clouds 
and its dynamics. In the disk of typical spiral galaxy, the magnetic field strength is usually 
estimated to be of several to more than 10/zG while in some regions of spira l galaxies the mag - 
netic field strength may be higher than several tens of microgauss (see B lBecket al. 1996h . 
By that estimations and others, it is almost certain that MHD density waves should also 
play an important role in the dynamic and evolution of structures within a magnetized gas 
disk (IFan fc Lou 19961 ). As is well known in classical nebular hypothesis by Kant-Laplace 
the condensation in protoplanetary rotating disk plays an important part in forming stars 
and planets. That part of the Kant-Laplace theory remains valid also for a contemporary 
cosmogony. The "standard" theory of t he multistage accretionary formation of planets, or 
the so-called core accretion mechanism (jSavronov 19721 : iPollack et al. 19961) r emained the 
most popular until recently, when it was criticized by (iBoss 20021 : iBoss 20031 ) and others. 
The main problem of the latter is t he timescale , which is longer than estimat es of the life- 
time of many planet-forming disks (ITaylor 19621 : iFeigelson fc Montmerle 19991 ) . In any case 
all theories rely on instabilities as a mechanism to transform a relatively uniform rotating 
gaseous disk into a planetary system. That is, at an early stage, the protosolar nebula are 
formed by fragments that separated from a molecular cloud. Planetary formation is thought 
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to start with inelastically colliding gaseous and dust particles settling to the central plane 
of the rotating nebula to form a thin layer around the plane. During the early evolution 
of the disk it is believed that the dust particles coagulate also into comets-planetesimals. 
On attaining a certain critical thickness (and, correspondingly, low temperature) small in 
comparison with the size of the disk, as a result of a local gravitational collapse the nebula 
disintegrats into the central body an d a number of separate protoplanets. Instab ilities arise 



as the thickness of the disk is reduced (IGurevich fc Chernin 19781 ; ISavronov 1972T ). If a rotat- 



ing gaseous disk has a large vertical thickness owing to a high int ernal temperature, then it i s 



sta bilzed again st gravitational instabilities by thermal motion (IGurevich fc Chernin 19781 ). 



In ([Boss 20041 ) it is demonstrated that convective cooling is able to cool the disk midplane 
at the desired rate to produce clumps in marginally unstable disks. The physical phenom- 
ena treated in the current paper occur during the stage of evolution of the protoplanetary 
cloud when the dust and gas in the disk start to condense into planetesimals and a star with 
current luminosity emerges at the center of the nebula. 

The rest of paper is organized as follows. In Sec. 2.1 we present basic equations and 
state our assumptions. In Sec. 2. 2 we present the dispersion relation to be solved. In Sec. 2. 3 
we pay particular attention to the conditions of the existing of complex-conjugate roots of 
the dispersion relation. We present our conclusions and discussion in Sec. 3. 



2. Hall MHD equations and the dispersion relation 

2.1. Basic equations 

We consider a thin rotating gaseous disk with angular velocities £l(r), where G is the 
gravitational constant and M is the mass of the central body, and r is the distance from the 
center of the rotating disk. The thickness of the disk can be estimated by c s /Q(r) where c s is 
the sound speed. The disk is made of partially ionized plasma where ions are well coupled to 
the neutrals while the electrons are not. However, charge neutrality, n e = rii, is assumed to 
be valid. The disk is immersed in a magnetic fiel d directed along the rotation axis (defined as 



the z-axis in our frame of reference) . Following (iBraginskii 19651 ) , the equations that govern 
the evolution of the two fluid system, namely the heavy particles (ions and atoms) and the 
electrons, are: 

the momentum equation 

dv -? 1^ -j GM _ 

nrrii— = —Vp H — J x B — nrrii — 5— r , (1) 
at c r 6 

here n is the number density of the heavy particles (ions and atoms), and the generalized 
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Ohm's law 

du -* e _ -> 1 - ^* ^ Vn e ej , . 

m e — = eE + -v x B jxB + T e - — , 2 

at c n e c n e an 

where <jr is the electrical conductivity. The generalized Ohm's law as given in Eq.(j2]) differs 
from the corresponding equation of MHD theory by the term on the left-hand side that 
describes the effect of electron inertia, by the third term on the right-hand side, which is 
the Hall effect, by the fourth term, which describes the effect of the electrons' pressure and 
by the last term that represents the drag force acting on the electrons. In addition, it is 
convenient to write the induction equation by substituting the electric field from Eq. (T5]) into 
Faraday's equation which than becomes 

— * — * — * 

dB -* -*. -* ,j x B . m e c -* du A -> 

where rj = c 2 /Ati<jr is magnetic diffusivity. The relative importance of the various terms 
in Eq.Q may be investigated by considering appropriate length and time scales. Thus, 
consider, for example, the second term on the right-hand side of Eq.([3]) (the Hall term). If 
we assume that V ~ 1/L, j ~ cB/AttL (displacement current is neglected), where L is a 
typical length-scale of the density inhomogeneity, then that term will be of the same order 
as the convective electric field if li ~ \/(L where 

c cJra. 



71: 



In other words, in order for the Hall term to be important the length-scale L should be of 
the same order of magnitude or less than the ions inertial length divided by i/C where £ is 
ionization degree, i.e., 

L < ii/VC. (4) 

Similarly, the ratio of the electron inertia term [the last term on right-hand side of Eq.([5])] 
to the convective term is given by: 

U) c 1 



f2 e Anen e L ^Airnmi 

Thus, the electron's inertial term is important for inhomogeneities which are characterized 
by L fa c/u pe = £ e . As far as the time-scale is concerned, it is seen from the relationship 
va/ 4 — VC^i ( V A = B / y/Aimrrii is the Alfven velocity) that for the Hall term to be important 
the frequency should be u > \/C^j while the electron inertia should be retained if uj > Q e . 
Note that in case of u » y/C^li the electrons drift in the wave's electric field, while the ions 
are immobile. Thus, if y/C^i << u> << Q e the second term on the right-hand side of Eq.(|3]), 
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namely the Hall term, is the leading term. This is the reason to term such approximation 
HMHD and waves in that regime Hall waves. It is seen from Eq. ([3]) that in such case 

C 2 Qi V H d 

CO PS 



CO 



pi 

,2 



L 2 L 



where vud = v\/(£liL is the phase velocity of the Hall waves in the presence of density 
gradients whose scale length is L. Such waves exist merely due to electron drift in the electric 
field of the Hall waves. It should be noted from the discussion above that the conditions 
for the Hall term to be significant are more easily satisfied as the degree of ionization is 
decreasing. Finally, it should be noted that for low enough densities the Ohmic dissipation 
term in E q.(3) (the last term on th e right hand side) is negligible in comparison with the 



Hal l term (jBalbus fc Terquem 20011 ). Additional support to that point of view is provided 
by ( Jin 19961 ) who estimated the ratio of the rotation period to Ohmic dissipation time to be 



of the order of \0~ 3 k 2 H 2 . As will be seen in the subsequent sections, the growth rates of the 
Hall instability are of the order of the rotation period and the relevant wave lengths satisfy 
kH < 1. Hence, the effect of Ohmic dissipation will be neglected from now on bearing in 
mind that it may nonetheless lower the growth rates of the investigated instabilities. 



2.2. The linearized equations and the dispersion relation 

We consider a differentially rotating disk with angular velocity Q(r) where r is the 
distance from the disk's center, and r-dependent density p(r). The disk is threaded by 
and axial magnetic field B = B (r)z. The equations governing the linear stability of the 
rotating disk may be derived from Eqs.(l)-(3) by assuming that the perturbations of the 
steady rotation are of the following form 

f(r, 0) = f(r)exp[i(m9 — cot)] (5) 

where f(r,0) stands for the perturbation of any of the physical variables that describe the 
system. It should be noted that strictly speaking expression (J5J) may be used only for rigid 
rotation, as differential rotation results in non exponential perturbations. Nonetheless that 
fact accentuates the nondependence of the Hall instability on the rotation shear, in contrast 
to MHD instabilities like MRI's. A full description of the temporal evolution of perturbations 
in differentially rotating disk will be described in a forthcoming publication. The amplitudes 
of the perturbed ^-component of the velocity, density and ^-component of the magnetic 
field can be expressed in terms of the amplitude of the radial component of the perturbed 
velocity. This results in an ordinary differential equation for u r (r) that should be solved with 
appropriate boundary conditions thus yielding an eigenvalue problem. However for purposes 
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of demonstration we first consider the simplified case of local approximation by assuming 

1 du. 
u r dr 



that k r <C k where k r ~ — ^ and k ~ — . After linearization, Eqs.(l)-(3) assume the 



following form: 

i(w - mD)ur + 2ttu e = (6) 
/ n . -X 2 rn 2 a m 2 6 Z 

( W - mfi)n e + z-u r = -c.- + -1^- (7) 

1 m Vr Bga + m^_ [w _ 55 i)] ^ L = (g) 



Cpfij r 87r p r Lr Bq 

(uj-mQ) w e = (9) 

p r 

where a is the perturbation of the density, and x is the epicyclic frequency given by x = 
+ 2rVLdVl/dr. The system of equations (6)- (9) yields following dispersion relation: 

- Co 2 u Hd - u[ X 2 + k 2 (c 2 s + V 2 A )} + u Hd ( X 2 + k2 ^) = ° ( 10 ) 

where V is the total unperturbed pressure, u = uj — mfi, k = m/r, and ujud = m ' j t^ L is the 
Hall drift frequency. It is a direct result of the assumed form of the perturbation, i.e., Eq.(5) 
that the axial and radial components of the perturbed magnetic field as well as the axial 
component of the perturbed velocity are zero. In addition, the derivatives of the equilibrium 
profiles have been neglected in deriving the linearized equations above, except in the axial 
component of Faraday's law (Eq.(8)) where, due to the Hall term, they are of the same order 
of the rest of the terms. 



2.3. The instability against non axisymmetric perturbations 

In the case of homogeneous density and magnetic field strength (L — > oo) the two roots 
of Eq. fllOp represent two stable branches of density waves that originate due to both the 
rotation of the disk as well as the external magnetic field. However in the case of density or 
magnetic field inhomogeneity the roots of Eq. ffTOj) with real coefficients are real if, and only 
if, the following conditions are satisfied: 

D = ^|(27X 2 + 4X(1 - 9F 2 ) - AY 2 + 8F 4 - 4F 6 ) < (11) 



where 



k 2 LV r V k 2 (c 2 + V 2 ) 



X = ^ + 
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Y 2 = x__ + k 2 (4 + v 2 A ) 



and 



. If the last condition is not fulfilled, Eq. fflO]) has two complex-conjugate roots one of which 
signifies instability. The conditions for that to happen are: 



I ! 2 ,,2 



X > — (-1 + 9Y 2 + VI + 9F 2 + 27F 4 + 27F 6 ) (12) 

— ( 



or 



X < —(-1 + 9Y 2 - VT+ 9F 2 + 27F 4 + 27F 6 ) (13) 

Further insight into the onset of the Hall instability may be gained by denoting the total 
pressure radial derivative by (p + Bq / 8ir) / Lj> where L-p is the inhomogeneity length of the 
total pressure. Next, the well known relation c s = HQ is recalled, and finally normalizing 
frequencies to Q, and wavelengths to H the following dispersion relation is obtained: 

Co 3 - qau 2 - Co[f + q\l + P~ 2 )} + qa[x 2 + £q\l + /T a '/2)\ = (14) 

where 



a 



(3 = c s /Va, q = kH, x — x/^-> an d ^ = L/L-p. It is first noted that for a — > the MHD 
regime is recovered and the roots of Eq. (IT41) represent the stable combination of the fast 
magnetosonic waves and the epicyclic oscillations. However, as a is increased (which means 
that L is decreased relative to the inertial length of the ions) the system enters into the Hall 
MHD regime. In that case, elementary analysis of Cardano's solution of cubic equations 
reveals that the nature of the roots of the dispersion equation (fill) hinges upon the value of 
\x which is given by: 

_ x 2 + eg 2 (i + r 2 /2) 
" x 2 + g 2 (i + /?" 2 ) ' 1 ] 

and the roots of the following quadratic equation: 

4/iS 2 + (1 + 18^-27^)5 + 4 = 0, (16) 

where 

q 2 a 2 



X 2 + g 2 (l + r 2 )' 

Thus, Eq. (fT4|) has two complex roots and hence the system is unstable in the following two 



cases: 
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1. \i > 1, for S\ < S < S%, where Si and S2 are the roots of Eq. ([TBI) . 

In this case £ must be positive which means that the density and the total pressure 
change radially in the same direction. It is therefore obvious that regions of instability 
occur where the total pressure changes more rapidly than the density (£ > 1). This 
is indeed the case in polytropic disks for which L/L p = 7 > 1 where L p is the inho- 
mogeneity length associated with the pressure. Hence £ > 7 and consequently // > 1, 
which means that the disk is unstable under the Hall instability if (3 is such that S is 
between the two roots of Eq. (TiT)]) . Exact values of the growth rate, obtained by the 
numerical solution of Eq. (fTlj) . are depicted in Fig. 1 for Keplerian rotation (x = 1), 
with l i jL\fC s = 10, and £ = 5/3, for various values of /3. It is found that the disk is 
unstable for 1 < (3 < 20. 

2. /i < 0, for S > max(Si, S 2 ). 

In this case it is obvious that the gradients of the density and the total pressure must 
have opposite signs (i.e. £ < 0). Such situations may occur when radial inflow plays an 
import ant role in the dyn amics of the evolving disk, such as in young protoplanetary 



clouds (jHogerheijde 20031 ). or when radial rings of non mono tonic density profiles are 



formed due to gravitational instabilities (IMayer et al. 20051 ). In these cases, in the 



limit a ^> 1 one of the solutions of Eq.( fT4l) is approximated by equating the second 
and fourth terms in Eq. fTH]) and is given by: 

u = ±i ll (17) 

where 

7 = ^\i\q\l + ^/2)-x\ (18) 

It is clear that in this case, the rotation plays a stabilizing role. Also, since S oc 1/(3, 
there is an upper bound on (3 for instability to occur but not a lower bound. Further 
more, in the limit of small (3 the growth rate grows without bound as (3 is decreased. 
Numerical solutions of Eq. (TT4|) for Keplerian rotation, and for iijL^JX = 10 and 
£ = —1.5 are depicted in Fig. 2. The instability exists for (3 < 9 and the growth rate 
is indeed a growing function of 1/(3. 

The unstable branch of the dispersion relation (TbTj) . termed the Hall instability, has been 



investigated in detail in (ILiverts fc Mond 20041 ) where it has been shown that it represents a 



quasi-electrostatic slow mode in which the perturbation in the ions' density and velocity play 
a crucial role while the perturbed electrons' density and magnetic field are negligible. Thus, 
the Hall instability provides a powerful mechanism for the azimuthal breaking of radially 
stratified disks into small fragments of size comparable to the disk's thickness and on time 
scale of the order of the rotation period. 
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3. Summary and discussion 

This paper examined the instability of a weakly ionized, thin disks threaded by an 
external magnetic field within the Hall-MHD model. The vertical stratification as well as 
the azimuthal variations of the disks properties were ignored whereas radial distribution was 
considered. In particular, in astrophysical context such weakly ionized gaseous nebulae are 
rele vant to the protoplanetary disks. The cond itions in protoplanetary disks were discussed 



in (jSavronov 19721 ; iGurevich fc Chernin 19781 ). At a certain stage of its evolution the star 



nebula is believed to have a characteristic disk size up to order of 30-100 AU and the total 
mass of the disk is believed to be less than roughly 0.1 of the mass of the central star. This 
yields an integrated column density of S ~ 3 x 10 2 g ■ cm~ 2 . Assuming that the mean mass of 
the particles is m p = Am#, the number density is given by n = E/Am#if (r) ~ 2 x 10 14 cm~ 3 
where H(r) is the scale thickness of the disk. To estimate the value of the factor A it should 
be noted that within this system, lighter elements such as hydrogen and helium were driven 
out of the central regions by star wind and radiation pressure during a highly active phase, 
leaving behind heavier elements like Na, Al, and K and dust particles. Thus in the outer 
part of the star nebula, ice and volatile gases were able to survive. As a result, the inner 
planets are formed of minerals, while the outer planets are more gaseous or icy. Concerning 
thickness estimation one can use c s (r)/Q(r), which yields H(r) pa 0.002AU. It should be 
noted also that due to very low temperature of the protoplanetary disks the only sources of 
ionization are non-thermal, e.g., cosmic rays, X-rays and the decay of radioactive elements. 



Following (ISano &: Stone 2002bl ) the ionization degree at the midplane of gaseous disk is 
estimated as £ = n e /n pa 10 -12 . So the disk material is a partially ionized plasma where 
ions and charged small dust grains are well coupled to the neutrals but electrons are not. 
Following this assumption we can estimate that tij \/C, is up to order of 1 AU. It should be 
noted that such big values are obtained mostly due to low degree of ionization, however one 
should keep in mind that existence of positivly charged grain particles increases the inertial 
length (see definition after Eq.(3)) and thus enhances the effect of the Hall term. In Sec. 2 
it has been demonstrated that the Hall term is important if H(r) < L < lijyfC,. Thus, due 
to the very small values of the ionization degree in protoplanetary disks the HMHD model 
has to be employed when studying the stability of structures with realistic radial density 
inhomogeneities. Indeed, radial stratification with length scale L > H may exist in the disk 
due to such mechanisms as axisymmetric density waves that give rise to alternating high and 
l ow density rings. On the o ther hand, L is bounded from below by H due to thermal pressure 



( IGurevich fc Chernin 19781 ). Hence, as has been further shown in Sec. 2, protoplanetary disks 
with such radial density distributions are susceptible to strong nonaxisymmetric instabilities 
whose growth rates are of the order of the rotation period of the disk. Such instabilities result 
in breaking of the density rings into fragments that may be identified as planetesimals. 
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Following widely adopted standpoint, an accumulation of planetesimals may lead to the 
next stage of evolution of protoplanetary disk, which is the coalesce of the planetesimals into 
protoplanets. Such planetesimals may survive t he thermal pressure if the ir characteristic 



size, i.e., 1/k, is bigger the disk's thickness H (IGurevich Sz Chernin 19781 ). On the other 
hand, the linear analysis presented in Sec. 2. 2 is valid if kL < 1. Combining those two 
conditions, and taking into account condition (j3J) results in the following limitations on k : 

< kH < 1 (19) 

where value of f3 has been taken as 1 for simplicity. Thus, for typical values in protoplanetary 
disks, namely, C = 10" 12 , ^ = 10 4 x B Q /Gs~ l , Q m 1.9 x lO'V 1 the left hand side of the 
second inequality in ( fl9l) is of the order of or smaller than unity. It is therefore again the 
small ionization degree that enablles the onset of the Hall instability, in the small magnetic 
field limit, and by thus providing a mechanism of initiation of the standart scenario of planet 
formation. 
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Fig. 1. — Growth rates of the Hall instability for Keplerain rotation, and £i/Ly/( = 10 and 
£ = 1.66. 



